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A Singularity Free Analytical Solution Of 
Artificial Satellite Motion With Drag 

by 

G. Soheifelo, A. Mueller and S. Starke 

GENERAL 


This report is broken down into three parts. 

PART I (by G. Scheifele) gives the connection between 
the existing Delaunay-Similar (DS) and Poincare-Similar (PS) 
satellite theories in the true anomaly version for the Jg 
perturbation and the new drag approach. It also gives an 
overall description of the concept of the approach. The ne- 
cessary expansions and the procedure to arrive at the computer 
program for the canonical forces is then outlined in detail. 

PART II (by A. Mueller) describes the procedure for the 
analytical integration of the equations developed in PART I, 
In addition, some numerical results are given. 

— PA.RT III (by S. Starke) describes and documents the com- 


puter program for the algebraic multiplication of the fourier 
series which creates the FORTRAN coding in an automatic manner. 
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PART I 

THE EQUATIONS OF MOTION FOR THE DRAG PROBLEM 
AND THEIR EXPANSION INTO FOURIER SERIES 

by 

G. Scheifele 


1. INTRODUCTION 

The objective of the theory described in this report is 
to arrive at a fully analytical theory for the motion of an 
artificial satellite which is perturbed by the term of 

the zonal geopotential expansion and by a drag force which is 
tangential to the orbit and proportional to the square of the 
velocity magnitude. The density function 

C(Xi ,X2 ,X3 , t) (I.l) 

does not have to be specified for developing the theory at 
hand. This is probably one of the most significant features 
of the theory. It is achieved by postponing the process of 
creating the final computer source - language coding, which 
is FORTRAN - compilable, to the very last stage. In this last 
stage (step No. 2 below) this coding for the final expressions 
is created in an automatic manner by a minicomputer. The 
corresponding program is written in BASIC language which al- 
lows easy alphanumeric string manipulation. This program may 
also be executed on a large computer if a BASIC compiler is 
available (as for instance on the UNIVAC 1110). 

The two-body elements used are of Poincare-type (PS- 
elements in tie true anomaly). These elements have been de- 
drived from the DS (Delaunay-Similar) elements which were 
first presented in Reference 1. The corresponding regular set 
of elements (PS-elements) has been derived in Reference 2, 
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and the corresponding differential equations for the perturbed 
motion are given in Reference 3. 

The drag theory is built on top of the J ^-theory in PS 
elements which is the subject of References 9 and 10. Only 
the first order short period and secular perturbation of the 
Jg-theory have been taken into account. The long period ef- 
fects of Jg and of the higher harmonics have been omitted 
because they are of the same order of magnitude as the uncer- 
tainiuies in the drag force and density model. Some coupling 
of drag and Jg is obtained implicitly because the total 
energy is used as a canonical variable, and is evaluated for 
the drag forces. 

In general, the theory described here is very similar to 
the one carried out by Lane (Reference 7) and Lane and Cranford 
(Reference 8) except for the fact that those approaches are 
based on the classical Brouwer- -orl approach (Reference 6) which 
uses classical elements (time as independent variable), while 
here the new PS formalism of perturbed two-body motion is 
used (true anomaly as the independent variable). In addition, 
the density model used here can be different from the ones used 
by Lane and by Brouwer-Hori and is not a fixed input to the theory. 

No reference is made to noncanonical approaches to the 
solution of the drag problem. The author considers these non- 
canonical approaches to be not adequate because they do not 
make use of the powerful tools which are typically provided 
by all hamilton mechanics approaches to orbital mechanics 
problems . 

Canonical treatment of the drag effects is done by using 
the transformation rules for the "canonical forces" as outlined 
in References 4 and 5. In order to allow an analytical inte- 
gration of these canonical forces they must be expanded into 
fourier series with respect to the true longitude which 

is one of the canonical angle variables in the PS-theory. 

These expansions are carried out in two steps : 
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l&t step. Manual computation of the fourler expansions 
of individual expressions which are not too complicated to 
expand. These are for instance the positive and negative 
powers of r (as they may occur in (I.l)), expansions of expres- 
sions of the type r^v, r^v^ and of the derivatives of the 
cartesian coordinates and of the time with respect to the PS 
elements. These expansions are carried out in such a way 
that terms of order magnitude OCe**) (e=eccentricity ) are 
maintained. Subsequently each of these expansions has been 
tested out on a minicomputer to verify the expected convergence 
properties . 

2nd Step. Or a WANG 2200 minicomputer a program was estab- 
lished which automatically discards higher order terms, carries 
out the multiplication logic of a product of up to four fourler 
series, and then produces a FORTRAN compilable alphanumeric out- 
put. This "custom tailored" algebraic processor is described 
and documented in PART III of this report. 

This two step procedure has the advantages that, first the 
errors are avoided which are created by manually multiplying 
fourier series, secondly the theory is flexible in the sense 
that step number two can be reproduced at any later time to 
either include a different model for the density and drag 
force or to increase or reduce the number of terms accounted 
for in the expansions which were carried out in the first step. 

The resulting theory is of relatively concise form and 
the resulting FORTRAN program compiles on the UNIVAC 1110 
EXEC 8 system and it is executable in the interactive mode. 

Some numeviaal results are given in PART II, where the 
analytical integration procedure is outlined explicitly. 

We can conclude that the use of the new Poincare-Similar 
elements in the true anomaly version is very well suited for 
developing a concise near earth satellite theory that accounts 
for a suitable density model. An attempt will be made to 
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further reduce the coding of the resulting FORTRAN program 
and to incorporate some accurate density models. 

2 , THE BASIC SET OF PS ELEMENTS USED 

The true anomaly DS elements are transformed into the 
singularity free PS elements by the following canonical 
transformation 


(moment as) 
a = p = 

s 1 

0 , = p = /2(<I'-Gycos(g+h) 
6 2 

0_ = p = /2(G-H)‘cosh 

^ 3 

''b = P* = ^ 


(coordinates) 

0^ = ()>+g+h 

0 = - y2($-G)sin(g+h) 

' /-7 0 

0 =„v/2(G-H)'sinh 

3 

0 = A 


The DS-elements are interpreted as follows; 

G; total angular momentum 

H; X -component of angular momentum 

3 

L: total energy 

related to two-body energy 

canonically conjugated elements 

g: argument of perigee 

h: argument of ascending node 

i: time element 

<t>: canonical true anomaly 

With q being 


q 



2/217 




2/^ 


(1.3) 


the transformation from the time t to the new independent 
variable t reads 

dt 


dT 


q 


(1.4) 



-13- 


The hamiltonian for the perturbed two-body mtjtlon is then 


F 


Pi- 



(1.5) 


(V=perturbing potential) 

and the initial conditions must be chosen such that F will 
initially vanish. 

From Reference 2 we can record the following abbreviations 
and expressions which will be used later. 


G - 

» = 1 

■(P 3 +P 3 ) 

(I. 

. 6 ) 

G = 

Pi- 

|(fj| + P 2 ) 

(I. 

.7) 

H 

G 

cosl 


(I. 

. 8 ) 

e== = 

II 

(— *^ - p(p| + P2) 

(P 2 +P 2 ) (1 

.9) 



J 





= Pa cosoi - Oj sinPi 

Z = pz sinoi + Or cosci 
2 


( 1 . 10 ) 
( 1 . 11 ) 

( 1 . 12 ) 

(1.13) 

(1.14) 

(1.15) 



-14- 


e sin(}> >= Q Zz 


(1,16) 


o cos<|> « Q Zi 


R = a, cosai + p 3 sinai 


( 1 . 17 ) 

(1.18) 


By using the above relations the transformations from car- 
tesian coordinates and velocity components to PS variables 
can be written as follows: 


xi = - ^ cosoi 


X2 = - + r sinoi 


X 3 


= -^/l-^(a| + pl)‘ H = || /2 Tg+hT 


/ t : 4 g 


" • 1 G ^ 3 , 

Xi=r (cosoi - - r(— 2 - sinOi + R) 

• • 1 G P 3 . 

X2 = r (sinoi - -^paR) + cosoi - 2^ R) 


^3 = (r R + r R) -|?i /l -;^(a| + p!)' = (rR)‘ /2(G+H)‘ 

/G 'io 2 g 


(1.19) 

( 1 . 20 ) 
( 1 . 21 ) 

(1.23) 


(1.24) 


where r, r and R are given by 


- = 1 + Q (p 2 COsai “ 02 sinai) 


O 1 * ^ 

3 . = (V ~ + G) (a2COSOi + Pasinoi) 


R = (pscosoj - Oasinoi) 


(1.25) 

(1.26) 
(1.27) 


The last transformation equation that reeds to be given 
here is the one for the transformation of the physical time t 

t = Oi, + ( 2 x!) 3 /a (^~‘t>~ ^/ 1 -e^'e sin([)) (1.28) 



with 


E-4' - - urctK 


sin4* 






(1.29) 


('f)sv 


unci whori' c siiitji , t' coH<f. :iri* 1 iikon 1 rom (1.15), 

(1.16), and U.17). 


3 , TRANSFQ RMA TIQN OF THE C ANONIC AL FOR CES 

Tht* forniul itsiii nl t.r:tu.sl()ru^ln^^ additional i'oroo.s which 
c’.an not bi^ dorivi-d Irom a pol.i?ntiul I'tinction under canonical 
variable transiormat. ion.; hue i’irst bin»n glvt?n by Brouwer and 
Horl (Roieronce 6). It has boon |j,-.jneral izi»d to inc'dude canon 
ical ayKtems which arc based f>n indt^pende-at variables difl'er- 
ent i'rom titm* in Reicrcncc a. A thoroufth and exhaustive de- 
scription oi this ex(. elision can be lound in Ri^t'erence 5 
(Chapter VIII). 

Let us tio back tcj tin.* cviir.u'al I'ortnul at ion ol the per- 
turbed two-body problem ii> . a ‘ v*.sian counii nates where the 
canonical equations of irn’'tiou read 





i)F‘' 



(It 

‘'I'k 







dt 


is the 

h amiltoiiiau 

ot the 

per 

extended 

phase space: 




- ,7(p-i 

+ P? 

P3 


X. 




(k-^1 .2,3,4) 


(1.30) 


(1.31) 


a 

r 


(1.32) 


where u is the gravitational parameter, V is the perturb- 
ing potential of Uic part of tlu^ aeopotenl ial expansion 

and pi( is the negative total energy: 
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V 

L 


|(p! + p| + p1) - + V 


“1 


e -- 
^ 3 


^3 2 1 

<-) -i 


3. J p 
2 *^2 <© 


(1.33) 

(1.34) 

(1.35) 


tR^: me^n radius of earth). 

The canonical forces of the x-type, Xo , Xi , Xz, X 3 , 
are absent at this initial stage and the canonical forces 
Pi. Pe. P 3 are the original forces given by the drag: 


Pj = C(x^^,t)v Vj , (d-1,2,3) (1.36) 

C(xj^,t) is the product of the density with the ballistic num- 
ber of the body. In this report, C is restricted to be a 
function of r cnly, namely the Laurent series given in (I.l). 

According to the rules given in References 4 and 5, a 
canonical force P 4 must be introduced which is defined by 


= Jl 


X - ^ P ) 
i i aPi 


(1.37) 


The left hand part of this expression is zero and for the 
right hand part we obtain 


P 4 - - piPi - P 2 P 2 - P 3 P 3 - C v^ (1.38) 


(observe that (1.30) and (1.31) imply Pj^ = Vj^) . 

At this point we are ready to transform the independent 
variable from t to t as posted in (1.4). Again, accord- 
ing to References 4 and 5 we can proceed as follows: 



(1.39) 
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-0 , (k>=l,2,3,4). 


(1.40) 


In the next atep we will i<ivo the transformation of the 
canonical lorces under the canonical transformation which 
ioads from the cartesian variables to the PS-elements . The 
transiox'mation .sclKiiiie can b.^ described by 


^k 



->• p 

•\t 


’\4 


\ 

' 'i’k • 


> u 


(1.41) 


and the explicit formulas are 


IJP], 

'^1 " k^l ^ 


(.1 = 1, 2, 3,4) 


^ 9Pu 'V 
^ kSi ^ 


(1.42) 


(1.43) 


They reduce to 


1, .X, c»X 

S P — - 
k=l ^k 


(, r\,, 3X, 


“j - kEi ■’k 


(1.44) 


and by inserting tlie expressions (.1.40) and (1.36) for P^ 
we obtain 


Tj = CM (V ,£,v^ 


“3 “ t CCO (V kS, 

(Observe the relation x,^ = t) . 


3x^ 

8t 

- V ^ 


8p . 
1 


(j=i,: 


9t 

- V^ 

80 . 

J 

30. 

2 


(1.45) 


(1.46) 
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ThesG equations form the basis for all of the subsequent theory 
which will deal with the explicit expressions and with the ex- 
pansions of and . 

The final form of the equations of motion is 


dT 


8F 



(k=l,2,3,4) 


dT 


3F 


+ U, 
k 


(1.47) 


For the case ~ = 0 (no drag) t ’ equations (1,47) 

have been solved. This is the J^-perturbed solution 
of the motion of an artificial satellite described in Refer- 
ences 10 and 9. The corresponding equations have been, pro- 
grammed and tested out (Program. PSANS). This Jg-solution 
is not discussed here and we will proceed to the expansion 
of the canonical forces Ti , Ta , T 3 , T 4 , Ui, Uz, U 3 , U 4 into 
fourier series. 

4 . FOURIER EXPANSION OF THE CANONICAL FORCES 

To enable an integration by quadrature of the canonical 
system (1.47), we have to expand the canonical foces and 

into fourier series with respect to the true longitude 
The prr^cess to arrive at this expansion is rather lengthy, and 
we will only outline the major steps and give the results. No 
intermediate calculations will be given here. 

Throughout the expansions, terms of order magnitude 0(e‘') 
will be maintained. In some cases where derivatives of the 
expanded expressions will be needed, the terms had to be ex- 
panded to order magnitude 0 (e®). 
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Let us list some additional abbreviations for quantities 
which occur very often: 


s^ 

= |(oi 

+ pl) 

« O(e^) 


(1.48) 



> 

(2P4)^'^“ 

3 = 


(1.49) 


y 


/i? 




= 3Zi = 

S(P2. 

cosci - Pz sinPi ) 

= 0(e) 

(1.50) 

Cz 

= 6Zz = 

3(P2 

sinoi + Pz cosp 1 ) 

= 0(e) 

(1.61) 


= s® 0® 

= 0(e = 

^) 


(1.53) 


The quantities Ci, C 2 and n are dimensionless (their 
physical dimension is of unity) . 

For the velocity magnitude v the expression 

Pi. = - I v2 + H - V , (1.53) 

which stems from F=0 is used, giving 

v"" = 2(-p., + ^ “ V) (1.54) 

The expression for r in PS variables is 


V and r , as \/ell as products of mixed powers of v and 
r will be expanded later. 

Let us first concentrate on the terms 


3 3xj^ 

k^l 9p^ 


V.h 


V, 


9x^ 

yo” 


(j=l,2,3,4) (1.56) 
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which occur in (1.45) and (1.46). After having inserted 
the PS elements into the expressions (1.56) and by making 
extensive use of the formulas (1.19) to (1.27), a consid- 
erable number of cancellations will occur. The result is an 
astonishingly simple final result for the above eight expres- 
sions. The explicit algebra to arrive at these simplified 
expressions takes several pages of hand-computation even if 
presented in a compressed form. These derivations are not 
presented here, we should like, however, to point out that 
the resulting formulas have been tested for correctness with 
a computer program that compares the original expressions 
with the ones obtained after the lengthy calculations. By 
taking into account the results of these algebraic manipula- 
tions we arrive at a new form of the canonical equations, 
which will serve as a starting basis for the fourier expansions 


T. = ^C(x)< 


0 

0 

+ V 

QZ 

^(1-3^6=") 

2 q 


2 ^ 
_ 0 _ 

3Pj 



(j=l,2,3,4) 



(1.57) 


U, 




G 






0 

1 „ 

^ 3 a. 


- V ^ 

at 1 

3 o. { 

[ 

2 

J 




0 






The expressions which need to be expanded into fourier 
series are now 


r^ V, r^ v^ , C(x ) , 


3r Sr at at 

aa.' 3p.' 3a.' ap, 

2 2 1 j 


(1,58) 


We will now proceed to these expansions. Expanding r and 
its powers is straightforward. As an example of these 
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binomial expansions we will give the one for r : 
r = p(l-ZiQ + - Z?Q3 + ZiQ" - Z?Q=) + 0(e®) (1.59) 


The expansions for Q are found by binomial expansion too, 
resulting in 


Q 

= 3 / 1 - 

s^6^ 

2 “ 

3(1 

. - 

1 

2 - ^ 0(e^ 

(1.60) 


= 

I 

) 


(exact) 

(1.61) 

Q' 

= 3^(1- 

I 

+ 

3 

32 

s"3”) 

+ 0(e®) 

(1.62) 

Q** 

= 3*‘(l - 

s^3^ + 

1 

4 

a** 

3") 

(exact) 

(1.63) 

Q' 

= 3^(1 - 

■ I 

4* 

15 

32 

s'^p-) 

+ 0(e®) 

(1.64) 


The fourier expansions of Z^i however, are becoming rather 
complicated since Z‘i , for instance, will create twelve terms 
which are powers of O 2 , Pz , sinci , ccsci , which each have 
to be expanded into individual fourier polynomials. For this 
reason it was decided to expand all the expressions in (1.58) 
into power series with the general term 

3^ P2 oz Zi zi sinoi cosci , (1.65) 

and to multiply these power series at the end in an automatic 
way with each other, establishing a table for the fourier 
series of expressif-ms of the type (1.65). By observing the 
relations 


Z? 

= 2s^ - 

Zi 

(I 

.66) 

Zi 

= 2s=^Zi 

- ZiZi 

(I 

.67) 

ISI 

= 4s“ - 

4s=^Z| + Z^ , 

(I. 

.68) 
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the highest power of Zz which will occur can be restricted 
to j«l. 

If these general guidelines are followed, the* expressions 
for the quantities in (1.58) can be expanded with respect to 
the eccentricity to obtain power series with the general term 
(1,65). Where these expansions are more or less straightfor- 
ward no intei’mediate steps are given. Here again, the formula 
manipulation to arrive at these expressions is generally very 
tedious and lengthy. For this reason, the resulting expres- 
sions have been programmed and their expected convergence 
properties have been tested out numerically by comparing with 
the exact expressions (1.65). The dimensionless quantities 
and ?2 are used. 

Expression for r^v: 

r^v - /2pi,' j^l + ^ t 15 n^) + ^f(2 + 27 n^) 

(1.69) 

+ I - X ^ ^ ?i)] + 0(e = ) + 0(V.e=^) 


Some remarks are in order concerning the last term in this 
expression. By evaluating the drag on a J 2 orbit the jorder 
of magnitude of the potential V will be of order magnitude 

_3 

Jg, i.e. 10 . Since we are only taking into account terms 

of order magnitude 0(6“^), and we are limiting the eccen- 
tricity to the interval 0<e<0.1, terms which .iave V as a 
common factor need only to be expanded to order magnitude 
0(e), thus giving four significant digits in the limit case 
e=0.1. The expression for V can be obtained from the con- 
dition that the hamiltonian must vanish in extended phase 
space ; 



= 1 + 2Zi Q - p? 8“* - 2Zi Q p? 0** . (1.70) 
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Because this expression contains information beyond the third 
digit only, it can be simplified by using 

Q = 3 /l - I = 3 + 0(e=^), (1.71) 

i.e. replacing Q by 3 , resulting in 


V 


P't 


= 1 - p? 3" + 2Ci(l - p! 3“) + 0(e") 


2 o4 


(1.72) 


By observing in addition 


= 3^y 


(1.73) 


and neglecting O(e^) terms we obtain the final expression: 


r®v = I y 3^ J 


4j2(l + n^')® - 1 + p? 3'i] 


- 2 Ci(2 + 15 + 2 Pi3‘^) + 2 C?(2 + 27 n^) (1.74) 


+ 4 d - 17 Cl 


+ O(e^) 


Expression for r^v^; 

Similarily, we obtain for r^v^: 


r^v3 =-|3^y f);. 


4 p? B“(l - p? S'* + 8 n'' + 8 n**) 


+ 2 p! 3'‘(2 + 2 p! 3" - n^) 

+ 2 e?[^pf 3“(-4 + 6 n^^) - 2 - 15 
+ 4 d(l + 2 p! 3") 


(1.75) 


- 9 P? 3" U O(e^) 
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3r 


ExpTeaaionB for and 


Sr 

Sa 


i 


The expressions for the derivatives of r are found by 
evaluating equation (1.59). This equation was given to the 
fifth order accuracy because taking the partial derivatives 
will reduce the order by one. Using the expanded versions 
of the Q - powers, which are given in (1,60) to (1,64) we 

3 T* 

can write a general expression for — : 




Zi 3(i“|n=^) + z? 3^(i-|n^) 


z? 4. 


- p 


- p 


3Q_ 


SZ 


Zi - 


2Z? 


sa-|n' 


) + 3Z? 3^ - 


4Z!f 




2Zi 




(1.76) 


+ 3Zf 3^(l-|n^) - 4Zj 3“* + 5Z^ 3j + 0(e ) 

(j=l,2,3,4) 

9 1!* 

and an equivalent expression holds for 

The eight partial derivatives of 3 (they are needed 
for the Q - derivatives), of Q, of p and of Zi which 
occur in the above expression are listed below. Derivatives 
of these expressions which are not listed are vanishing. 


93 ,, 3 

9p4 4Pn 

|2j = - Oa 


(1.77) 

(1.78) 

(1.79) 


9Q 

9p4 


_3_ 

4L 


(1,80) 
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3P - 

2p 

^ 

3 P z 

U3^ 

lIL- . 

2a 2 

3 a 2 


l£_ = 

V 

3p*t 

2p5 

3Zi 



(1-n*) 


= - z- 


8a : 


(1,81 


(1.82) 


(1.83) 


(I .84) 


3Zi 

= oosoi (1,85) 

9pa 


3Zi 

= -sinai (1.86) 

902 


Inserting these expressions into (1,76) yields, after the 
truncation of higher order terms and collection of the terms 
has been carried out, the following expressions: 


9r 

3pi 


= 0 , 


Sr 
3p • 


= 0 , 


= 0 

3a, 


Sr 

3ai 


= 0 


(1,87) 


3r 
3 p 2 


2p 2 




- Ci(l-§n^) + 5?(l-|n^) - r,? + ?'{ 


+ P 023'[t.i(|+ - I i;? + I d “ c'i] 


- p Scoso 


n 1 


(1.88) 


^ n" - 25i(l-|n') + 3c!a-|ri=') 


n - 


- 4 + 5 cl 

1 


+ O(e^) . 
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||-= - [l . + 5*] 

(I. 89) 

“ 4?T ■ 2t;ln-y]^) + 3ci - + 0(e®) 

If^ = p Sz [i “ I - -^ o'* - 2?i(i-|n=^) + 3cf(i-|n=^) 

(1.90) 

- 4?f + 50^ + 0(e®) 


s!(l-fn^) - ?! ?'] 

+ p p 26 “[ci(|+ in“) - I ?" + I c! - ?!] 

(1.91) 

+ p BsinCTi |l - § ^ n** - 2Ci(l-.|n^) +3?i(l- 

- 4d + scij + 0(e®) 

3 1 Bt 

Expressions for -r— and — ; 

30k 3Pk 

Let us first record the time equation (1.28) again, and 
write it by using the new abbreviations. For this purpose 
we may first list some basic relations 

/l-e^' = 1 - (a|+p|) = 1 - , (1.92) 

2y 


e sinoi = Q Za 


, e cosCi = Q Zi 


(1.93) 


ifilco 
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Thua the time equation may ruiw be written in tlu* form 


t 


Oi, + 


2p.,3 = 


CE-4>) 


r 

P 


(1-tr) Q 


Lj Z;. 


where atanda for 

ZzQ 

E _ (j, « _ 2 arctan 

2“T)^+ZiQ 


(1.94) 


(1.95) 


Let us introduce the quantity T 




(I.9G) 


and first expand E - ^ , yinee this is the difference be- 
tween the true and the ecetuitric anomaly it is of the order 
of the eccentricity, This implies that; the art>;ument of the 
arctan function is small, allowinj? us to expand it into a 
power series about the point zero. The result Is 


E - (j) = - Z?3(l-iriM I’ + ^ 


+ 0(c’^) 


(1.97) 


In the above and all the subsequc-nt expansions for t we 
have to be (’areful to take intti account terms up to and in- 
cluding 0(e'') , because tif tlu’ derivatives to be taken later. 
Inserting (1.97) into the time equation (1,94) yields 


t 


a 


+ 


1 

2p 



r ‘4- 


1 

12 


r ^ 
1 2 


A. 

80 


~(l~n') ZzQ 


(1.98) 


+ 0(eM 
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The last term —(l“n^)Z 2 Q is expanded by inserting ^ from 
(1.59) and then using equations (1.60) to (1.64) for the 
powers of Q . The result is, after adequate truncations to 
the order needed : 


^-(i-nMZaQ = ?2 1 - ?i(i-|n'‘) + d(i-Tn') 


- d + Cl - I + I n”] + 0 (eM 


(1.99) 


We may now write down a general expression for the derivatives 

St 

which equally applies to the derivatives : 


j 


j 


3t 


904 


3li 


9p^ (2p4)^" 


- Kz r(i-|n") + I’' 


- - I - Ci(i- |n^) + c? - 


Cz [1 - I - 2 Ci(l-|r|2j + 3 ^2 „ 4 


(2L) 




9Ca 


r 5 2.1 

[1 - 4 " -"4 


n" - Ci(i-|nM + Ci(i-^Ti^) - Cl + c"} 


(I .100) 


+ r(i-^n'')-i ^2 r'(i-|n") + ^ Cz r = 


+ >. 9(n^)/5 I 2 3-.7„2.1p l-2p3\ 

+ Cz C 4 - 2 h ^4 I 6 ^ 


- [1 - f - 1 


(j=l,2,3,4) 
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In tho rurLher proet'duro tho poworn oT r ,2 urc* roducod to 
n“l by using tho relation 


rj * 2 n® - r/: 


( 1 . 101 ) 


which Is equivalent to (1.66). The following expressions 
need to be inserted in (1,100). The results of their ex- 
pressions are given below. 


r = 1 + I n* + I r- - I 5ia+fn^+||ri') 


+ I dU+nO - + i c! - 


r2~ X + ri^ + f n" 


- 2 Cl + Xg Cl 


Ci(i+ |n^) + 




p3 ^ 


1 + |(n^+n“) - Ci(|+-^n^) + Ci(|+3ti®) 


3 j. 13 - 4 

4 ^ 16 


r"^ 1 + 2 n'’ + I n“ - Ci(2+|n") + cfcl + ^n^') 


5 , 3 j, 35 _ ii 

- 2 " 10 


(1.102) 


(1.103) 


(1.104) 


(1.105) 


cSr^= 4 n" - 4 c! r/ + g 


2 K. 2 J. " 


(1.106) 



-so- 


ar 




3p . 

J 


+ 0 ?; (| + ^n^) - I c! + 


j|^!] 


(1.107) 


1 9(n^) n. /3 . 15 2 \ lf-2.|. 5 a 1 a"! 

ap, >- -J 


Last, the derivatives of ?i,C 2 and need to be listed. 

The ones which vanish are omitted in the list. 


3?i 

9pi 

9Ci 

ap4 

3? 2 

302 

3 (11=" ) 

3p2 


" ?2 ) 


Cl 

4p4 


3 cosa 


P2 S' 


3Ci 

3Ci 

= 3 COS0I 


9P2 

30 2 


9C2 

3 C 2 

11 

1 

1 

301 

3p2 



3C2 

1 ) 

— ' — 


3p4 

3 (nM 

3( n^) 

= 02 3' , 

302 

3pi* 


= - 3 coso'i (1.108) 


= 3 sinoi (1.109) 


(I. 110) 

4p<4 


(I. Ill) 

2p„ 


At this stage, all the information is provided to insert the 
necessary quantities into (I. 100). However, at this stage too, 
the individual expansions for the r" and for the derivatives 
have been tested out numerically before proceeding with the 
algebra. After a considerable amount of explicit computations 
the following final expansions for the partial derivatives of 
the time with respect to the PS variables result. 
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6 n^(l-ri^) + Ci(4-13n=') - C?(6-21n^) 


+8 Cl - 10 C** I + 0(e®) 


3t r 

JC 2 sinoi 12 (1 - n^)- ?i(20-29n^) 

Scfj IGpi, 3 I I— 

+ 2S Cl - 36 Ci^ 

+ cosai|l 6 - 12 n^+ 2 ri "-4 Ci( 3 - 4 ri^) ( 1 . 113 ) 

+ c!(12-19ti^) - 12 Cl +12 cf] 

- C2 02 3 Js - 5 - 10 Cl + 12 CiHI + O(e^) 


(1.114) 


(1.115) 


0 


(1.116) 
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at 11 
a p 2 16 p I* 3 

- sinoi 


?2 cosaijla (1-n^) - Ci(20-29n=^) -t- 28 „ 36 ^3^ 




12 + 2 n** - 4 ?i(3-4n^) 


+ (12 - 19 n*") - 12 cj + 12 ?i 


0 


(1.117) 


+ ? P 3 

Z 2 


8-5r| “10?^+12e 


2 

1_ 


at 

9p3 


+ 0(e®) 


= 0 


(1.118) 


at 


1 C: 


a p 4 16 p 4 3 


io2 


20 - 7 n® - Si (12-611=^) + 8 C? - 4 ci 


+ 0(e=) (1.119) 

Now all the necessary derivatives and expressions are 
expanded which will be needed to evaluate (1.57) except for 
the density function C , To test out the theory and the 
general procedure outline in the next section we will insert 
a constant density for the preliminary numerical tests. Var- 
ious expansions of simple representations of the density 
function are currently being investigated. The general form 
of the resulting expression for a spherically symmetric den- 
sity function is an expression which is very similar to the 
ones for the partial derivatives of r , and only powers of 
will occur. If latitude and time dependent density func- 
tions are considered, the resulting expressions become more 
complicated, and will be dependent on powers of ?2 « 

5 . AUTOMATIC GENERATION OF THE FINAL EQUATIONS 

As it was outlined in the introduction, this theory is 
of a "dynamic" character. This allows for example a continuous 
updating of the density model. It is of course, necessary to 
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develop such a function into an expansion in powers of Ci> 
52 , sinCi, cosCi by explio-it manual computation. Also the 
expansion must be complete to a certain order of the eccen- 
tricity; in our case 0(e‘*) terms must be included. Then, 
before the expansion may be used for an input to the auto- 
matic algebraic multiplier its expected numerical convergence 
behavior must be tested out. 

Let us here record the expressions (1,57) for the canon- 
ical forces 
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which can be simplified to 

X 

u = (4-5n‘^) + 0(e®) (1.124) 

4p2p 


needs to be considered. 

left 

By denoting the left hand terms by , the middle 

terms by and the right hand terms by 

IC K. 

write the following multiplication table 


rnXef t 

k 


(r="v).(C)- 



(k=l,2,3,4,5,6,7,8) (1,125) 


T 


T 


mid 

2 

mid 

4 

mid 

5 

mid 


(r'^v) ' (|~) • (u) *(C) 


(r=^v)‘(|^)-(u)-(C) 


(r=‘v).(|~)-(u).(C) 


(r^v).(||^).(u).(C) 


^right ^ (j.2v^-(|~)*(C) 


, right ^ 


= (r^^v^) •(!—) ‘(C) 


(1.126) 

(1.127) 

(1.128) 

(1 . 129) 

(1.130) 

(1.131) 


mright 

5 


(1.132) 
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right ^ (1,2^3 
u ^ ^ 2 


(1.133) 


^right ^ (r2v®).(C) 


(1.134) 


The terms which are not listed here are absent, as can 
be seen by inspection of the expressions developed in the last 
section . 

It is of interest to note, that due to(I.125), and 

are equal up to a different common factor; is the 

only canonical force which gets contributions from all three 
terms in (1.123). 

Of special interest is the force Tj^. It is the only 
vanishing one. The consequence is that the variable oi 
(true longitude) is not affected by the drag perturbation , 

This is a peculiarity of the PS-theory. Even though this 
seems strange on a first look, it must be emphasized that it 
does by no means imply that there would be no in-track effect. 
This in-track or timing effect of the moving body within the 
orbit is accounted for by T^ , which determines the motion 
of the time element Oi, . This interesting observation re- 
flects the fact, that the geometry of the motion is fully 
separated from the dynamics within the orbit, which is typ- 
ical for the new canonical variable formulation of the DS 
and PS theories in the true or eccentric anomaly case. 

The ten products displayed in (1.125) to (1.134) are 
carried out with the automatic algebraic multiplier which is 
documented in PART III of this report. The task to be per- 
formed by this multiplier is twofold: 

(1) It will determine if a term is of order magnitude 
0(e®) or higher. If this is the case, the term 
is discarded. Practical experience show that 
about 85% of the terms occuring in these products 
will be discarded. 
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(2) If a non-discardable term is encountered, it is 
identified, and a pointer is established .ihich 
refers to one of the elements listed in tnu 
"Table” below (Tl, T2, •••, T27) . The perti- 
nent factors as well as the pointer are then 
temporarily stored and saved on tape in such a 
manner that the resulting statements are FORTRAN 
compilable. 

The »TABLE": 


A careful inspection of the expressions occuring in 
(1.125) to (1.134) and of the expansions of their product 
terms which were given in Section 3 shows that the only terms 
of order magnitude 0(e'*) and smaller are the ones given in 
the following "TABLE". The expressions for the right hand 
sides are computed in a separate FORTRAN subroutine and are 
transferred by means of a common block. 
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TABLE 


order 1 . ignitude 0(1) terms: 

T22 = 1 
T26 = sinoi 
T27 - cooai 

order magnitude 0(e) terms: 

= T1 = 3C-0Z sinOx + pz oosOx) 


Cl 

sinOx 

= T2 = 3(- 

1 

2 

Pz ^ Pz sin^Ox 



+ — 
2 

Pz 

aos2ox ) 

C 1 

aosOx 

= T3 = 3(| 

Pz 

1 . o 

- Pz s^n 2 pi 



+ ^ 
^ 2 

Pz 

G032Pi ) 

Cz 


= T23 = 3(pz 

sinOx + Pz ecePi ) 

Cz 

sina 1 

= T24 = 8(| 

Pz 

1 

+ p Pz sun2ox 



1 

2 

Pz 

ooe2ai ) 

Cz 

0080 X 

= T25 = 6(| 

Pz 

+ 5 Pz sin2oi 




Pz 

oos2ox ) 


(1.135) 

(1.136) 

(1.137) 

(1.138) 

(1.139) 

(1.140) 


order magnitude 
Cl Cz = T4 

Z.X ?2 sinOx = T5 

Cl Cz eosoi = T6 

cl = T7 


O(e^) terms: 


= 3^ U (p2-Oz) si«2ai 
+ pz Oz eos2a 

Pi l1 1 

= 3^ - f sinax + ^(Pz-o|)oosai 

+ ^ pzCfz sin3o\ - i( pz-ol )cos3o i^ 

= 3^jj (Pz-Pz) sinOx + | PzPz cospi 

+ ^ (Pz-Oz )ain3a 1 + pzOzCosSo,^ 

~ 2 (Pz~^z) ~ pzOz sinZox 

+ ^ (pz-Pz) cos2aiJ 


(1.141) 

(1.142) 

(1.143) 

(1.144) 
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Si sinoi 


= T 8 =02 


1 1 
j <p2+3ai) oinoi - paOz eoedi 


1 H 

4 (Pz-^l) ein3a p p cos3a 

E l ^ 1 ^ ^ 

• p2ff2 sinai + -^(Spl+cri) coppi 


(1,145) 


(1.146) 


1 1 
- 2 P 2 CI 2 sinScTi + -^(Pz-al) £3os3a^ 


order magnitude 0(e ) terms: 




= TIO = 3^ T P2(p|+cf2) sinOi + ~ c^2(p2+cr|) cosCi 
*- (1.147) 

+ 7 Pz (pI+Sp?) sin 3 a +7 a_(3p2-p^) cps3p~l 

1 


Cl ?2 sinoi = Til = 0^ ^ PzCPz+cjf) -^02 (pl-crt) sin2ai 


1 ** 2 2 2 

.2 ,,2 


1_ 


1 *1 

- 2 Pzcrf pps2pi + ^ p^(3p2-a|) stn4Pi 


(1.148) 


1 

R p 2 (pl-30^) cos4p 
a z l_ 


Cl C 2 coscTi = T12 = 0®j^ a 2 (p|+o|) + i p 2 (p|-a|) sin2pi 

1 X (1.149) 

+ 2 *^2a2 eos2pi + g P2(pl“3p2) sin^Oi 


+ <Jz(3p|-a2) eos4ai 


r 3 

Cl 


= T13 ” 3^ ~ § Cf 2 (p 2 +cr|) sinoi + “ P 2 (pz+ci 2 ) ctp.'icfi 


4 

1 

4 


aa ( 3p2-pi , sin3o\ 


(1.150) 


-t 

+ -T p2(P2-3ai) aos3a 


i_ 


Cl s'unai 


= T14 = 0^ 


“ 3 

J 8 


OziPz+Oz) + Pz(Pz+3a|) ein2oi 


+ ^ of aos 20 i + P2(pz-3p|) sin4pi 


(1.151) 


+ aa(3p|-a|) cos4p^^ 


n 


Cl oosoi 


= T15 = 0® 


O -1 

g P2 (Pz+< 7|) - ^ az(3p|+p|) sin20i 


1 , X , , (1.152) 

+ 75 - p| pos 2 pi - -Q a2(3p2-02) sin4ai 


1 1 
+ g Pz(pl-3a|) cps4ai 


L) 
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order magnitude 0(e'*) terms: 




= T16 = (Pa-cJa) ain2ai + — PzOzipl+al) cob2o\ 

+ |j (P 2 +O 2 ) - I PjfJj (1.153) 


+ ^ PzOzip'i-Oz) oob40i 


Zl Ca Binai = T17 = 3 " 




p o (p|+a|) ainai + (P 2 -a 2 ) cosOi 

22 o 


+ '2 P 2 P 2 sin3oi - 


(p^-3an 


_ (1.154) 

3 1 1 

+ 3 p|ai eos3oi + ^ P 2 O 2 (p|- 0 i) sinSCi 
-[^ (P 2 +O 2 ) - I pIpIJcosSoiJ 
52 eosoj = T18 = 0'*|§ (P 2 -o^) sinOi + i P 2 '? 2 (pl+P 2 ) aoaai 

+ j^ (3p'i-a|) - I pla^cinSoi 
+ I pIpz 003301 +1^ (p^+a^) 

- -| pfc^sinSai + ^ pzOzCpi-al) t?()c5pij 


(1.155) 




= T19 = 3 " 


I (pI+cjI)^ - P202(pi+a?) sin20i 


1 I (1.156) 

+ 2 tp2“^2) ^^ 0320 ^ - 2 p202 (P2-02) sin 4 o i 


(P 2 +P 2 ) - T p 2 aiU’o. 34 ai 


5i 8^nOl 


= T20 = 3*' 


{ i^p 


5+5aa ) + 7 P 2 P 2 i ainoi 


^ t'o.34a i| 

i]si, 


- 2 P2CT2(P2+02) OOSOi + 


(3p^-5on 


3 1 1 (1.157) 

+ g p|a| sin3oi + ^ p2a2(P2+3a|) posSoi 


16 


(P 2 +CJ 2 ) - i Pfoll stnSai 


'0 

+ ” PaOa (p 2 -a|) <?os5aij 



Qosai 
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= T21 a 3'*|- i p2CJa(pi+o^) sinai +jj (SpHoU 

+ I pfo^oosai - I p2p£ (3pi{-ai) sin3ai 

+ [^ (5p^-3a|) - I p|a^oc7fl3CTi (1.158) 

- “ P2a2(p|-ff|) aindoi +|^ (p2+a^) 

- § pIp|1cos5pi| 
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PART II 

ANALYTICAL INTEGRATION PROCEDURE 
by 

A. Mueller 


1 . INTRODUCTION 

The subsequent considerations are concerned with the 
analytical integration of the equations of motion defined in 
PART I . These equations have the form 



dT 



dx 


3F 



i-1,2,3,4 



T 


i+A 


(II. 1) 


where P is the hamiltonian which includes the two-body and 
Jg potential and the drag perturbation is included in the 
canonical forces T^ . Because of the coupling considerations 
given in PART I we will assume that the integration may be 
separated into two integrals. 




' 

T^ dx 


dp. 

— i dx = - 

dP 

dx + 

dx 

3a 


‘i+4 


dx 


(II. 2) 


The integration of the part due to the hamiltonian has 
already been treated in Reference 9 and we proceed to the 
integration of the canonical forces. 
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2 . CANONICAL FORCES 

As described in PART I, the canonical forces have 

been expanded about the eccentricity so that these forces may 
be written as a fourier series in the fast variable , 
which in turn will facilitate the integration. However, to 
analytically integrate the forces, they must be expressed as 
a function of the independent variable and the initial con- 
ditions. As is standard in perturbation theory, this is done 
by inserting the two-body expressions for the motion of the 
elements into the equations for the canonical forces, These 
two-body expressions are obtained from the integration of the 
canonical differential equations due to the two-body part of 
the hamiltonian. It is known that this integration is a 
oanoniaal tvanefovmation in itself, thus to be consistent, 
the canonical forces must undergo another transformation. 

The transformation in the elements is defined by the 
following equations 


ai(T) 

om(t) 


T + Ci(0) 


u 


(2p,) 


3^ 


T + Cl, (0) 


(II. 3) 


c^(t) = for k=2,3 


P^(t) = P^CO) for k=l,2,3,4 


As defined in Reference 5 the canonical transformation 
of the forces is given by 


T. (t , a ) = S 

j ’ i k=o 


9p, 


3c. 


T, 


3p^(T) 


+ T 


k+4 




(II .4a) 
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and the conjugate forces 




= , r 

k=o 


3P, 


3a, 


+ T, 


^ 3a, (T) 3a, (T) 

j j 


(II. 4b) 


considering (11,3) one finds that the transformation reduces 
to the following 


Tj(T,ajJ = for j=l,2.3,5,6,7,8 

(11,5) 

T 4 (T,a, ) = T., + -^■ ^ • —7 - T Te 

(2pOV^ 

The canonical forces are now in a form to be easily integrated. 
3 . INTEGRALS 

Let us define the right hand integral of (II. 2) as 


RiCT.aj^j 


T^(t,Cj^) dT 

J 


(II .6) 


Due to the fact that the canonical forces have been ex- 
pressed in a fourier series, the indefinite integral may 

be expressed by a summation of terms which are multiples of 
the following integrals 
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Porms marked by the arrow appear only in the integral Ri, 
and are a result oi the independent variable appearing in 
(11.5) as factor. 

PERTURBATIONS 

Partial coupling effects may be introduced by evaluating 
the integral.^ with the elements which are predicted 

from the Jg satellite theory. The procedure by which the 
integrals are evaluated and .he perturbations due to the 
and drag are summed requires four steps 


1. 

given 

T=0 , 

a^(o) ^ (II. 6) R^(o) 

2 . 

given 


(o) 

3. 

given 


. T -V (II. 6) R^(t) 

5 . 

given 

\(T) 

and update • 




= - [Rk(0 - Bk(o)] 



Pk CO 



The algorithm requires no iteration process to determine 
the value of the true anomaly as a function of both drag and 
the Jn perturbations as in Reference 8. This is due to the fact 
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that the anomaly is a canonical element and is only perturbed 
by the ob.' ateness pocential. 

5 . THE i^EAN M0TI0^^ 

In the integration of (II. 6) we considered the energy pj, 
to be fixed. However the drag force dissipates energy and 
thus the mean motion is changing with the independent vari- 
able. In the expression for in (1 1. 5) the error in 

the mean .motion is multiplied by t again, so that the resul- 
ting error for considering the mean motion fixed is of order 
Y being the relative magnitude of the drag force. 

This error may not be neglected for large t 

Let us assume that the energy dissipated by drag is a 
linear function of the independent variable. For large values 
of T this dissipated energy AL can be expressed as 


AL 


Rs(t) - Ra(o) 

T “ 

AL 

T 

L 

_At_ 


This is true because the .secular term dominates the periodic 
terms fer large t in Ro(t) . 


Also if the dissipated energy AL is small compared to 
the total energy L then one may express the mean motion 
term in (11,5) as 


3pT 

3)i 

-r 4. 5 

T T — 

— • — t 

AL 


(2(L-AL))'''' 

(2L)®'"" 

2L 

At 




- 48 - 


This new term introduces integrals of the form 
T ^ dx = T ® / 3 





2nx 




These terms are premultiplied by factors of order 7^ , Also 
it should be noted that the procedure to evaluate the drag 
perturbations on the solution remains unchanged except for the 
fact that perturbation of the energy el'^ment must be known 
before the perturbation of the time element may be computed. 
Also the additional terms due to the mean motion in the per- 
turbation of the time element do not require additional FORTRAN 
coding to be developed by the algebraic manipulator. From 
(II. 5 ) one can see that the expression for Tg may be used 
to compute the additional terms. 

The procedure outlined in this section corresponds to 
what is commonly known as "the second integration of the mean 
motion" . 
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6. NUMERICAL EXAMPLES 

To demonstrate the accuracy of the expansions and the 
Integration procedure, the analytical solution has been com- 
pared to a numerically integrated solution for several test 
orbits. Both the numerical and analytical techniques tssume 
a drag model of the form 


where C, = 2.2 

d. 

I = 100 Ibs/ft^ 
p = .5xl0~^ kg/m^ 

The density p , is to be a constant. The value con- 
sidered is the approximate density of the earth's atmosphere 

at the very low height of 175 km. The value of the ballistic 
M 

number, , is the average for the Shuttle. Also, both numer- 
ical and analytical models neglect the inertial velocity of 
the earth's atmosphere. 

Both the numerical and analytical solutions will include 
the perturbation due to the oblateness of the earth. 

Three different orbit cases (Table 1) with eccentricities 
of e=0, e=.015, e=.l, were chosen for integration. In Table 2 
the position difference between the numerical an>. analytical 
solutions are compared after 20 revolutions. Also displayed 
is the position difference if drag is neglected in the analyti- 
cal solution. 

In all three cases the analytical theory accounts very 
well for the large perturbation due to drag. The large posi- 
tion differences when drag is neglected is from the "in track" 
error. The results in Table 2 confirm that the basic analyti- 
cal formulation is sound. 
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a(km) 

e 

I 

w 

0 

M 

Case 1 

6678 

0.0 

0 

0 

0 

20° 

Case 2 

6678 

0.15 

30® 

0 

0 

20° 

Case 3 

7300 

0.1 



30° 

0 

0 

20° 


TABLB 1 

Orbit Cases (initial conditions) 


Model 

(D 

11 

O 

e=. 015 

e=. 1 

Position difference 

Neglect 

Drag 

1481 km 

1506 km 

1920 km 


.97 km 

1.01 km 

2.18 km 


TABLE 2 

Differences of analytical vs numerical integration 
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PART III 

THE A .GEBRAIC MULTIPLIER 
by 

S . Starke 

1 , INTRODUCTION 

The algebraic multiplier^ was programmed in BASIC on a 
Wang 2200-T minicomputer (16k bytes). Extensive use was made 
of the alphanumeric character string manipulation offered by 
BASIC. 

The program documented below has the capability of multi- 
plying up to four fourier series with respect to and each 

individual fourier series may be multiplied by a common factor 
AK$. The trigonometric terms are stored in a separate alpha- 
numeric array ZK$(J) , and the corresponding factors are 
given in array elements EK$(J) . 

The program -input consists of fourier series, which are 
stored on tape, and up to ten fourier series can be loaded 
into core sequentially to allow automatic (buffered) opera- 
tion of the multiplier without manual interaction. 

During execution products of the type 
Z1$(I)*Z2$(J)*Z3$(K)*Z4$(L) 

are searched for terms which are not of higher order. If a 
low order term is encountered, the program wilT identify the 
term by pointing to the ''Table" which is stored in memory. 
Subsequently, ASCII Coding for a FORTRAN compatible state- 
ment will be created and saved on an output tape. With a 
separate telecommunication program, the output tapes are 
dumped to the UNIVAC 1110 computer where they are compiled 
after some minor modifications have been made (Subroutine 
headings, dimension statements, common-bloacks ) . 


+ 

Multiplier and manipulator will be used interchangeably. 
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For the following description of the program the follow- 
ing reference should be consulted. 

"System 2200 A/B Reference Manual, 

Wang Laboratories, Inc,, 1974". 
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2 , PROGRAM DESCRIPTION 

The A variables and arrays (Analysis) 


A1 

= 

Al(l) 

+ 

A2(l) 

+ 

A3(l) 

+ 

A4(l) 


power 

of 

the 


term 

A2 


Al(2) 

+ 

A2(2) 

+ 

A3(2) 

+ 

A4(2) 

= 

power 

of 

the 

^2 

term 

A3 

= 

Al(3) 


A2(3) 

+ 

A3(3) 

+ 

A4(3) 

= 

power 

of 

the 

P2 

term 

A4 

= 

Al(4) 

+ 

A2(4) 

+ 

A3(4) 

+ 

A4(4) 


power 

of 

the 

^2 

term 

A5 

Not used 












A6 

= 

Al(6) 

•f 

A2(6) 


A3(6) 

+ 

A4(6) 

= 

power 

of 

sin 

'"l 

term 

A7 


Al(7) 

+ 

A2(7) 

+ 

A3(7) 

+ 

A4(7) 


power 

of 

cos 


term 


In all cases the Al(5)->-A4(5) terms are the power of the indi- 
vidual term being multiplied. Therefore, if these terms are 
added we have the power of the resulting product. This is 
usually the primary test for excluding a term. 

C, Cl, C2, C3, C4 

C1-J-C4 are the numerical constants associated with each 
of the four terms with C being the product of C1->C4. 

The general form of the term being analyzed is 

C p2 <^2 ^1 

this term is held in the Tl$( )-+-T4$( ) arrays, 

T$(27)13 this array holds pointers if a term is present after 
a particular multiplication. An”*" indicates that a 
term is present. The first two characters of this 
array tell ifaZ 2 or a ios is present. Distinc- 
tion on the power of Z^ is made by the entry point 
and exit point of the table. 

The general form of the arrays being manipulated is 
A1${T1$(1)E1$(1) + T1$(2)E1$(2) + . . . + T1$(10)E1$( 10) } = FACTOR 1 

A2${T2$(1)E2$(1) } 


A4${ 


} = FACTOR 4 



-56- 


Al$ contains a character string of up to 20 characters. 

t They are never changed Just carried along. 

A4$ 

Tl$(n) contains the terms being multiplied, 

* 

T4$(n) n=l,2,***,10 

El$(n) these terms are considered constants for the 
; Tl$(n)->-T4$(n) arrays. They may have up to 20 char- 

; acters . They are never manipulated only carried 

E4$(n) along. 

N(1)*''N(4) this array contains the number of terms in the 
Tl$( ) T4$C ) and El$( ) E4$( ) arrays. 

I,J,K,L these are the loop counters for the nested Do-Loop 
(FOR : NEXT Loops) used by the multiplier. Once in 
the multiplier section these values must never be 
changed, 

M Loop counter used within the multiplier. This value 

may be reused as necessary. 

R1,R2,R3,R4 

these are the bookkeeping variables for determining 
how many terms were examined, higher order etc. 

H;H(9);H$(9) :H1$(9,10) :H2$(9,10) 

this is the density buffer. It can contain 9 terms 
of the type being multiplied. 

H: buffer pointer (H=l->9) 

H(H) : array holding the (N(4) value) 

H$(H) : array holding the A4$ values. 

H1$(H,M); array holding the T4$(M) values. (M=1->H(H)^10) 

H2$(H,M): array holding the E4$(M) values. (M-l->-H(H)£lO) 

As can be seen, the density buffer will only load into the 
factor 4 array. 
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B;B(3) ;B1$(3,10) ;B2$(3,10) :B$(3) ;B1 

This is the factor buffer. It will hold up to 3 
factors of the type being multiplied. 

B buffer pointer when loading buffer 

B(B1) array holding the N(2) values 

array holding the T2$(M) values (M^l-^BCBl)^. 3) 

B2$(B1,M) array holding the E2$(M) values (M=1+B(B1)£ 3) 

B$(B1) array holding the A2$ value 
B1 buffer pointer when unloading buffer (Bl=l')-3) 

As can be seen this buffer can only be loaded 

into the factor 2 array 

Ml( ) this array holds the extra term generated if a 

power of Z 2 is found. 

M9 a flag if an extra term was generated. M9=0 if no 

term was generated, otherwise M9 points to the 
extra term stored in Ml( ) 

All the other variables are less important as they are used 
for general bookkeeping. 

SUBROUTINES (as they appear in the listing) 

'1 to '4 general print routines used to print out 
factors 1 to factor 4 

'37 subroutine to print the input variables 

'16 used to load the table into memory 

'17 to '20 loads, respectively, factor 1 to factor 4 into 
memory 

'22 loads the density buffer from tape 

'23 loads the factor buffer from tape 

'30 initial entry into the multiplier (the starting 

point ) 

'10 re-entry into the multiplier if it has stopped 

because of an error (usually an end of tape error). 
This is the manual re-entry point. 

'll loads/resets the density buffer into factor 4 
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'12 loads/resets the factor buffer into factor 2 

'35(G$,Z) used to analyze the Tl$( )->-T4$( ) arrays. The 

position being analyzed is given by Z ; the string 
being analyzed is given in G$. If a value is 
present its power is returned in A . 

'36 analyzes the G$ string for the presence of a 

sin or cos. If a sin is found A is set to 1, if 
a cos is found B is set to 1. 

'41 subroutine to make the substitution for a power of 

Z 2 • The extra term is temporarily stored in the 
Ml array . 

'42 subroutine to read back the extra term stored in 

the Ml array by '41. 

'43(F0) routine to create a FORTRAN COMPATIBLE statement 

from the input data. It works in conjunction with 
'44. The resulting statement is in FS$ and F9$. 

'44 routine to pack the Al$->-A4$ and El$( )-+E4$( ) 

arrays into F8$ and F9$ inserting an when 

necessary and deleting all unnecessary blanks. 

The storing of data onto tape is also done by '43 and '44. 

The Fl$8 string contains the character used for generating 

the continuation cards. When a "Z" is encountered it is 

reset to an @(in HEX a 1 added to @ = A). 

General Starting Procedures 

1. load program "MULTFF" 

2. load Density buffer (S.F. key 22) 

3. load factor buffer (S.F. key 23) 

4. load factor 1 (S.F. key 1) 

5. load factor 3 (S.F. key 3) 

6. press special function key "START” (S.F. key 30) 

7. when told insert "TABLE" tape and press the continue 
button on the console not the S.F. key marked "CONTINUE" 
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8. when prompted check paper, insert a free data tape and 
press the continue button on the console not the S.F. 
key marked "CONTINUE” 

9. if the program has an end of tape error (Error code 49) 

a) remove the full tape 

b) insert a fresh tape 

c) press S.F, key "CONTINUE" (S.F. key 10) 

In general the machine will issue a message if the opera- 
tor must supply some type of information. 
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start outer 
loop 

1=1 to N(4) 




Multiplier 
Phase 2 


Analysis of 
important terms 


C 


C ZiZ^pzazslnCTicoso 


Analysis of 
4th term T4$(l) 

Analysis^A4( ) 




DEFFN'35 

DEFFN'36 


Start 1st 
Inner Loop 
J=1 to N(3) 


Analysis of 
3rd term T3$(J) 

Analysis“>-A3( ) 


1 ; 



Start 2nd 
Inner Loop 
K=1 to N(2) 




Analysis of 
2nd term T2$(k) 

Analysis^A2( ) 


DEFFN ' 35 
DEFFN ' 36 




1 
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580 



JL Note: 610 starts the basic 

comparison of all data 

610) stored in Al( )->-A4( ) 
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DEFFN'42 



[Multiplier! 

[Phase 5 J 

Incrementation of 
Loops and Buffers 
















STOP 
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66 . 



DEFFN'43 

(FO) 



iConvcrt C and M 
to characters 
and place in 
string f8$ 


Subroutine creates a FORTRAN 
card from the given input 

FO : a flag telling if this is 

the 1st card of tlie multiplication 
= 1 create a card using A1->A4 
= 0 create a card using 
El$( )-^E4$( ) 


rC=ClAC2*C3*C4 •, 
M=Table location 




DEFFN'44 inserts an 
between the character strings 
and omits all unnecessary 
blanks, it also checks to 
see if an overflow of the 
F8 string will occur. If an 
overflow occurs the information 
is stored to tape, a continuation 
card is generated and the process 
continues . 


Return 
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VARIABLES USED IN ALGEBRAIC MULTIPLIER 

SYSTEM 2200 VARIABLE CHECK-OFF LIST 

PROGRAM NAME DATE 

VERSION PROGRAMMER 

SYSTEM 




A 

B 

c 

D 

E 

F 

Q 

H 1 

j 

K 

L 

M 

N 

0 

p 

0 

n 

1 

B 

u 

V 

w 

X 

V 

z 


. 









: 

1 








■ 

■ 






■ 

1 













x: 




■ 







1 


i 

2 

ii 










n 

1 





i 




1 

■ 


■ 


1 

3 










■ 

■ 





■ 

■ 

■ 

■ 


i 

■ 


■ 



<1 

i 

■ 

■ 

■ 

■ 






B 

1 

■ 

i 

■ 

1 

B 

■ 

■ 

m 

■ 

■ 


■ 



& 

i 

■ 

1 

1 

i 




i 

j 


B 

■ 

B 

i 

i 

i 

■ 

i 

■ 

B 

■ 


■ 



G 


■ 

i 

B 

1 




■ 


B 

i 

■ 

■ 

1 

■ 

■ 

B 

B 

■ 

1 

B 


B 


1 

n 



S 

■ 

i 




■ 

1 

■ 

■ 

■ 

B 

i 

■ 

S 

■ 

■ 

■ 

B 

i 


B 


! 

1 

H 


s 

i 

B 

■ 



~1 

1 

1 

B 

■ 

■ 

i 

■ 

1 

1 

■ 

i 

B 

B 

1 


■ 



II 


B 

■ 

i 


n 



1 

i 


1 

■ 

1 

B 

i 

i 

■ 

i 

■ 

1 

1 

i 

■ 


s 

B 


B 

i 

B 




_J 

i 

■ 

B 

1 

B 

B 

1 

■ 

■ 

■ 

11 

■ 

B 

i 

iia 




NUMERIC SCALARS 
FORMAT = MN 


NUMERIC ARRAYS 
FORMAT = MN( 


ALPHA NUMERIC SCALARS 
FORMAT = MNS 


ALPHA NUMERIC ARRAYS 
FORMAT = MNS( 


NOTE: 

0 = NON COMMON 

1 = COMMON DEFINED BY THIS 

MODULE 

2= COMMON DEFINED BY PRE- 
VIOUS MODULE 
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Listing oi BASIC Program for Algebraic Multiplier 


■ nt-GEL^RftlC HULTIPUIER NRI'fTEN BV S. STFIPKE 
-.iO F.Eri..-+^>«^.. . VERS10f^=':i flDfif''T£L' FOR JSC BV Ci. 5. MFIRCH 

1.0 ,C a#i<= 

'-•0 D I M H ‘t B / Hi: '-i ;j 20 , HI -f- < 20/ H2'»: >' P.< 5 > 2 ; P=0 ' PE^i -t >LNS 1 T't' BUFPER - 

. BU DIM .B''Js':'/.Bt':2:'20.. BiiC:ii/J?:'2£C.B24 <Z/P:>2tHB!=0.r;!iiL -rHCTOP .BLinftB^ 
C'lfi 084 £4/ CiP't'iS 

?0 DIM T4‘:2?:-'12.. fll$20/ H2420.. fl3'420.- H442Q/ 2tU/ G420.. OI 420 ,, Fit .'.iMij 
ri D I M V9i.iO,. V84S : I N I T FP r« VP4 

7C- D'll'l fl420.. Bi-20/ Ml-: 4/ 3 :\. F8464/ FPilB/ P420/ P146/ TP4±1 
80 DIM Tii-''.li-:O20.. 724''.i0>20/ ‘f 24 <10 >20/ T44<i0>20/ M<4> 

_£<0... DIM £14<10>2/.£24<1 Q>2.. E24<10>2.. E44<10>2/ Ei-.S 


iOO D 1 M R1 < ?' / R2 < 7' > / fli- < ? > / fl4 < 7 > : SELECT PR I f Tf 005 't 64 > ; PR I fTf hLJ< < 02 > : Ef JD 
120 DEFPN '00;P=4;PRnir HLJ<<OD> :QOSUB ■'i;GOSUC •'2:G0SUD j'2 :Grj5iji:! '4‘ R[£TIJR 
H 

iI-0 DEFFfJ'l ... 

14G PR I NTIJS I NG 3 0.- 1/ fHil ; FOR I =1T0 N < 1 > : GOSUB ' 27 < 1,. ’’1 L < I / Ei.f < I > ; rjEXT I 

”150 DEFFM''2 " ' ^ 

160 PR I tn USING 20/ 2/ FI24;F0R 1=1 TO N<2>; GOSUB '27<2/ 7 24< I E24.ajL>j HEJ^T . 
I; PRINT " N<2> ="/ N<2>.. HEX<ODOft> * RETURN 

170 DEFFN''2 ... 

■ ISO PRINTIJSnJG 20/ 2/ R2-4 : FOR 1=1 TO fJ<2>; GOSUB -'I V':!*-/ T24< I >/ E24< I > > ; NEXT 

190 DEFFri'4 

200 PRlfTfUSINQ 20/ 4/ H44 ; FOR .1=1 TO. fb'.4> ;GOSUB '2V <4/ T44< I >/ E44< I > >.; NEXT 
I ; PR I NT " r J < 4 = " / N < 4 > / HEX < ODOR > 

210 RFJTURN , . . 

220 REM - PRINT UTILITV RDUTirmS- 

„,.220 U.N1.T.< 20 >_Fi-L2/. 1 N.PUT. . " FT LB _ t jR!‘ 1E ■■..._5.T R < F i- ^SEL ECT PRINT 015< 72> : 

PR I r-JTUS I NG 240/ STR < F 4 < 1 > .. 2/ S > 

221 STR<F4<1>/ 1/ i>=HEX<C2> :STR<F4<1>/ LEN<F4<i> >+l>=STR<V94/ 1/ 10-LEN<F4<1 
> >> ;STR<F4<1>/ 10/ i>=HEX<OD> : RETURN 

240 KILE r^F^ME = "m#####" . .... . . 

250 DEFFrJ-'27<fl/ H4,. B4> ; PR I NT USING 40.. R.. I.- R4.. R,. 1/ B4 . RETURN 

...270 .RSI'i .-LOflELTlEriORTl- 

2S0 DLFFN'i6:STi:ir- "TRBLE CRSSETTE ?":DriTR LOlTD 'nRBLE":FQR I=1T0 27 ; DflTR 

LORD Ti.<l> ; NEXT I ; RETURN .. .... 

290 STOP "POSITIOr^ TRPC USIfJG < DRTFl LORD 'NRME' 2"; RETURN 
200 DEFrN'17; PRINT "EOU i":QOSUD 290 : DlTTft LORD N<1> .• R14 ; FOR I-1.T0 !.K1>;D. 
hlR LORD T14<I>/ E14<I> NJEXT I:PRirJT "LORDED" : RETURN 
. .210 DLFFir.iO^FSlfilJ'mLiL'LUjOSUB . HPAJ-RTFl LORLl-N '1237^a2^l . FJl!E..Jj=im.ij.i::z:t2 12^ . 

RT fl LORD T 2-i- < 1 > .. E24’ < I > : r JEXT I : PR 1 f TV " LOf JDEL- ” : F:ET UF:N 

220 DEFFN '19: PRINT "EQU 2"; GOSUB 290;DRT R LOfiD rJ<2;>. R24 ; FOR 1 =iT0..N<2> ; D 
m fl LORD T24 < I > .. E24 < 1 > : NEXT I . PR 1 f JT " LORDED " . RET URN 

220 DLFFN"20;PF:It^T "EOU 4"; GOSUB 290 ; DRTR LORL- fi<4>/ R44 ; FOR I=1T0 tJ<4:-.D 
RTB LORD T 4i >' I > ■ E4 i •' 1'. : t iE! 7f I PP I NT " l-.iTRI .r- f - » PF 'I i 'f 'M 
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I d ir'i C”ian 


.^Z5 .1<LF.I-N '22. H==D-i='F2 
“■40 H=H‘i i:DRTR LORI' H-'H!''.. H4''H:' ; FO!:: I=i'fO ; DHTR LOnD H14 *'H. ! X. H2:t <H.- I 

,';NE.\'7 I.FF.Itn R. . 1 F H=5'TH&N .260 : DRTR LORD " + GOTO 240 

245 DCF FN' 22 ; B=0 : PRINT "FRCTDR BUI-FEP"; ;GOSUD 290 : ir^PUT "LOFlD N1 TH M-JL=i : I 

.. F W?-;2THEN 250;WS»2 

250 b=G+i:DriTR LORD ; FOR I-iTO E:*'B>:DRTR LORD I B2'4--'B.. I 

•> T PP 7 ^JT p. : IF ~i=:r r.nro i nori »-4-" nn7n -'sn 

'260 PRINT "BUrFEP LORDED" : RETUFTJ 
:-S0 REM -MULTIPLIEK;- 

400 DEFFN'20:F;EI'1 - S'l RRT -; SELECT PRIN7 0O5''B4r' : PRirn HEX-:03:> : GOSUB •'if:: 

-£’~^jj::i0.r4Ui- i ... :1 1 ■ . na sun-.'-L 2 

4iO STOP "PRPER ?".GOSUB 220; ST OP "CORRECT 7RPE ?" 


415 DflTR SR'v’E BT''H> F 4- >'>: SELECT PRINT 015‘'72T : GOSUB 'O0;GOTO 420 

. .420- XCFF n REM . RESTi7F:;T- RETURN- £l.EfiFl-_SELECX FlIiI7JT...r^l5 CTl::- ^-VSl =S'TE C. 

FT '■:i:-'. 2.. ST ; IF poS':Fi*:i:-'=Frri=o THErj 422 

- 4 F4 . 2, POS 't F 4 >' 1 > =F F T -1 

422 GOSUL: 220: PRINT " CONTINURTIOH OF ;";VS4:.i" >";STOP "CORRECT TRPE ?" 

425 DRTft SnVE OPEfJ "TCONTi" 

- 420 MR. Rl. R2., R2. F4.. R1, R2,.- RTU-R4=0 ; T^-J.-|i;^L...-M.-C:i-;.-C:2.. 034-04=1 

440 MRT Hi=2EP ; r^iRT H2=ZEP : MRT R3 =2ER : MFIT R4=2ER : 1 N I T >' 20 r- F81.. F9t- : Q4 = " ICO 

= Ri! 

450 SELECT PRIfTi' 005 < 04 i'' ; PRINT "DO NOT DISTURB - PROGRRri MORKIfJG - " ; SE 

-LECT PRINT 01.5<72> ;G0SIJB-^'42<:i^ 

470 ON P GOTO 5S0,. 540,. 510.. 4S0 

■480 FOR 1=1 TO NC4..' - - — ■ 

4?0 FOR M=i TO 4;Q0SUD '25r7 45 ': I 5+2=+=M;-' ; R4''M:'>=R ; R t''.5:-'=R4':5:-‘+fl : NEXT M ; CO 
—UVER.T -ST F:UT4.^4ia.-1.4-jUJlCT-44..CI..-i :• ! i'l 1 0 -, C 4 


500 GOSUB ■'2S':T44''I:0 ;R4>'t.>=R;R4':?r'=E 

-J310 FOR J=i-TO NC2> — - 

520 F OR M=i TCI 4 : GOSUB ' 25 T24 <1 J 7 .. 5+2^=M > ; R2 C M T =R : R2 ^ 5 !' =R2 •: 5 T +R : r JEl^a M : CO 

. -NVERT STR t.’ T3:T < J 1. FJUM 4 T24i: 4T T T-TO C2 — — - - 

520 GOSUB '264 T24' 4 H :< :< ; R3' 4 6 :• =R ; R2 4 7 > =B 

_540-FCi|'4-44=l-4LCi_W4ST 

550 F OR M=i TO 4 ; GOSUB ■' 25 4. T 2‘$ 4 }4 T . 5+2=4:M > ; R2 4 M 7 =R ; R2 4 5 =R2 4 5 > +R : r JE^-4T M 

—550 X'ON VERT S-TR4 T 244 K4-.. -1.- NUN 4-T24 CK T > >-- TO - C2 

570 GOSUB '264T244K> > . R24S:'=fl : R247>=B 

5S0 FOR M=1 TO 4 : GOSUCi ' 25 4 T 1-4 4 L .. 5+2'*M > ; R1 4 M T =R ; HI 4 5 ri =Ri 4 5 +R : NEXT M : GO 

4 - 7 > ~ B 


BOO CONVERT STR 4 T14 4 L lU 1.. NUM 4 T 14 4 L ? TO 01 
010 IF R145>+r;245> +024-5?' +R4‘%5-r:i-4 THEN -S30 
012 IF E14 4L) = "41. -Ri+Rl'+=B'+"+=4V THEN 615 ; GOTO 616 

. -615 IF-R145>+R2>:5:-'+R245C'+R44S:0-l -THEN-S90 

626 C = 131 C'2=+: C 3: •+: C 4 

_iS2'0_ai=RiXl l+R2-ill2i^3LC12-+R'l^ 
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original PAGE lfc» 

SZm R2=Rl''2;'»+H2>'2’'+f;l ''£;‘+R4i;2> POOR CjUAf/IiV 

IF R2-C2THLN 649 ; &OSUB ..-'41 . . 

640 H3-R1 1 > +R2 r 1' > +R2 ■' 1 !' +R4 ii' > 

_E50 fl4=Rii;4:'+R2<4>+nic4;>+-fi4‘:4:< 

660 R6=Ri t: 6 +R2 6 +R1 *' 6 :< +R4 6 ;• : I F R662THU;N 670 : 011 = " ' 6 IN-'”: Z4= " S " ; GOTO 9 
._LQ - ' *” " 


670 R7=fll ■: 7 :< +R2 ^ 7 > +R7 •' V '• +R4 7 ^ : I F-' R762THEr J EOGi ; Gii;= '• ■' COS ' " : Z4-= " C " : GOTO 9 

.. A ^ 

6S0 IF R6+R7-C2 THEN 690 : Gli:-" 'SI N-FCOS'' 24:=" 2" ; GOTO picn 

.690__K::r'=" " ; IF A2=G THLH 7.00 ±.. 10 = "J2"- 

700 IF R6=0 THEN 710 : STR'^R-Jv 2.- l.'‘ = "S" : GOTO 720 
OIF R7=0 THEN.-Z2 Q:S;TRi.R4... 2., 1 ^ = " C " 


720 R=R1 : 7=0 ON H GOTO 720,. TI O, 720 ; 2=2 : IF fi=4THEN 72 fJ : H=27/6 Z=0 
725 REH - - TRBLE SERF.’CH - - 

720 FOR Fi=6+H -5 TO 6*R-2':1F STR''T4 1.. 2;''=S7R-;R4-,. 1., 2>THr-N 740 ; rJEXT M ' R1 
T="TREL.E ERROR" ; Z$=" f" ; GO'f 0 .910 

7^0 Zi=STR‘:T4 -'M;:'.. 2 j 1> : GOSUO '42>'0:';1F Z4 = " + "7HCr^ 92fi ■ GOTH c^r1r1 




890 F:±=Ri+i ; Gcn o 920: REN — HIGH ORDER — 

. 9^00 R2=R2+i,G0T0 920 ; REM JJQ. CONS'L 

910 R1-R2+1 : PRINT " — ERROR — UNEKF’ECTED Gil.; " RT '•.i :PRItnUSItJC 

0.. L.' K.’ J.' I ; !1=99 . GDSUD . '.42' I'O r; 

920 R4=R4+1:IF l■19C=0 THCfJ 920 ; GOSUB '42; GOTO 621 

..920_rifiT._m=2ER.:Cl=l:rJEXT L:I F F'=1TH E N IGl l 


940 r-iftT H2=ZER : C2=l : f JE;4T K;IF P=2THEN 1011 

.>50 MH7 R2=ZER : C2=l : rJC NT ...T_ : IF FZLlZTHEN^lOii 

960 HRT R4=2ER . C4=i ; NCNT I 

...j970. REri -- MULTIPLICATION FINISHED -- PF: I N'L RESULTS -- 
980 IF K9=0 THEN 1011 

_ 9 . 9 . 0 ._GQ. 8 :U 1 ;l. 7 ^ , 


1011 STFPiFei,. 1>=F11 ; STR<F81.. 7:-‘ = " > " : DR'l fl SAVE FSl.. F 91 : DflTR SAVE E^^D 

_1012 PRirJTU£lNG .i450,. F81.t j_FR.lNTUS.irjQ 1460, F9T : SELECT PPirJT_ 015P72:'- 

1012 PRINT HL>:''.0D>.; " TOTRL TERMS = 'L R4> Rl/’" : NC=” "/rIV '* ; UN« 


: CONST 


L. rCT — II 


R4-Rl~r:2-e:2....HD;: < OD :? JLF 2 tF::F=P".JHEN.CODOR 7 


1020 GOSUO 1070; IF H:''=10THErJ 1040 ; PR I fFf USING 240.- 01 : PRINT :PPINTUSirJG 20 
.. 4 .. R4 : ppurr HEIN' ODOR 7 . GOTO 425 

1040 PRirjTUSirJG 240.. G!l : F'F'irJT ; GOSUO ■■2;G0SUB '11 ; PRirJTlJSIfJG 20. 4. RJ-1 F'PI 

„.NT. hejo;od0r:u.qotd 425 ' 

1050 REM END OF MULTIPLIER SECTION — 

..1060 DEFPrj 'll . H=0 : REl'j — — .DErJSIT-V LORD RfJD RESET 

1070 H=H+1:1F H:-=10THEN 1090 : fJ'M : FOR 11=110 H>'H.-' ; 'f 41':M>=Hil':H.. HI . E4 
CMO =H21 ':FLJ'i:u r JEia' Jl. «41=HiCHl_: REXLUtli 


^10.80 L'EFFrT'12 : 81=0 ; PEM FRCTOR LORD AND RESET 

*1090 81=81+1; IF 81>=W9+lTHLrJ 1110 .N':2>=B*;E1 :i ;FOR: M=1T0 ; T 21 ';M:'=BlSt; 

B1 M 7 ; E21 < M =821 C Bi, M ;• : f JENT M . R24 =B1 1 Bi j . GOSUD ' 4 ; PETURfJ 
4.110 SELECT PRINT 005 ■: 64 >; STOP ■ "PRODUCT FimSHLD. LORD fJEW EPi ! AND - 'ST 
RRT ■' - " : ErJD 

.■_11.2£1 .EiE£T21122i!<lGl..-J^l.Lfl^.'RI CrSTRCn i. i -■ ■ Tu:- p-y-n t mpo -i -t -;n ■ Q =n 

1120 IF ZC=?THLrJ il40:fl=R/2 

-.1140 RETURN . - - - ... . . ..._ . 

1150 DEFFrJ-'26':Gi:‘ ;fl,. B=0 ; IF STF;>rGl.. 14.. 2''0"SIfJ"T HCTJ il60:fl=l 

1160 -IF STRCGl.- 17. 2C*0"CC!S-" TUCN-1170:B=1 _ 

1170 RETURN 

_1.220 . DE.EOLS41 : FJEM P01..1FF- n, - . 7.2 Pnt io n ■ MQ ^riQ-! -1 



1240 hl-;M9 . i:''=2 + i. ; C=-C -72- 

. 12^'0 2 ;’'=Hi : Fl±=fll +2 

1260 Fi 2 .' ft 2 '’“ 2 ' : kE'f UP'M 

1270 DLFI^'N ^42 : F:t.F’i +^F;ERC' E’-.TRR F'RDDUCT-'-^ 
12S0 C=l'!ia'i?.. i ;< : Fll=t’il';i'1?.' 2 :' : R2==r’ii«:M?.. 2:1 
-. 12 Sa.l! 2 =l EE 7 Hr 2 U-- 


page 

OF POOR 


1200 DCFFfJ '42':F0:'' ; IF F0=1THEN li; 00 ; ETR>'F61.- 1>=F 11 ; CONVERT C ^0 STR‘'FS1.- 7 

tr. 4 t ■ " ■ 

STR''F8i^« M TO STR<F84-.* SiX* : F£:=LEN''F85 :''-^i 

1220 SELECT F'FlXNT 0Ci5*:e4:< ; GOSUD '44i:Eil-CL:'X:-G0SUC . ■'44XE21CK2' > ; G0SU3 '45< 
H~'l' '' J ' F' ‘ GOi'tJL' 44 •' E41 '* 1 .* 

"•■ n rr l ^J ■-»u^r:^T i 7n P. 1 -3. ■ r:nNVER.T_ FTH STP >'F '='l.. 1 0. ,1.:i..^!l^^2! 

jL“-4n cnNVERT J'f 0 STR’'F i<l,- il.> I?'.* <'4^!' ; CONVERT ITO STF:''FP4v 12^ llv ''#i' 

iv=i0 12 ^=zi;STR':f:si. i4:'=sTftCFjt':i.:-^^ i:- — . - - 

l^^o rONVERT eiTO ST R>:F91.- i:^- CONVERT HTD STR>'F?1.. 16.. IX- 

1370 F)L'<DC':F1Jv Oi;> : IF— £TRCF-li. 6. 1> CHEXCSEL'THEN XX:-'G;F11=" e." ;GOXO .12S 

?- '-•f 1 7 SELECI F'FX[ fJT-OXb '-72 ..GUSUE-144 CHliC>— STF;iEE;iv_7J =ST FX FEllv .3 

t:± QOSDCr -44.:R2i :- 
1 3' S 1 G 06-U B ■■ 4 4 .' R2 1 :■' ; G L'6"UB ' 4 4 ■: 4 1 :• 

iXI 2 ST F: ■: F?S.. LEt J •: F?1 +1 •: ” : 004 =F‘SX- ; C!S-f =F?4' . GOTO 13'?G 
1390 ON K9'-l GOTO 1292. 1394.. 13 96 - - - 

1392 h:9=K9 : 1 ; 0G1=F91 ; STR'XjOF,. 9. 9 ;.=STF:i:F'91.. 9.. B'> ; INIT.'20:''F8T.. F91 ; GOTO 13'99 


1394 K3-<=^:9-*-l : STR>:001, 29,. 22:'=ST'P>:F91.. 22:' ; INIT ■'20:'F94-.. F9T ; GOTO 1399 

139C K9=K9H ; STR.'.QGIj 51.. 14;'=STRCFS1j3".. 14:- ;STF:-;Q91v 1, S:>=STRCF£-J£j^2-l^.S:< . - 
i390SELECT F'RINT 015.‘'72:' : F'RINTUSING 1450.001.; ;PRINTUSINQ 1460.. 091 : ORTR S 

ftVE 0.01-. 091 .INI T •; 20 :■ FOl .. Faiv-OSl.. QD4c_J?3=Q 

13’99 RETURN 

_1400..J:‘E.FF1L-45 >: RTC;-:- i f ... J~!7gn-3THE U ...1 40::3,J41;53^03: 


1402 K1=LEN<R1> ; IF F1=1THEN 1440: IF F8+F1X64THLN 1420 

1410 .S.TR>:FS1.- Fs:. = " + ’' :STR-:FE1. F£::1:'=P1.:FS=LEN';F£'1:'+1.QDT0 1440 . - 

1420 O01=F91 . C!3^1=F91 : GOSUO '47 : PRINT HE'-0-:07> : STP''F91.. f , i:. = ''4.'' ; PS=7 : GOTO 

..1410 - 

1440 RETURN 

1460 'I################ 

1500DEFFN ' 47: SELECT PRINT O±5':72i:'.;PRIJ.j:n,lSING-1450..Q&la-jPF.;INTUSING 1460.. G! 
91 ; C'RTH SRVE QOl.. 094 : !NIT':20;' F94.. F94-.. OOl.. 091 : K9=0 : RETURfJ 

-1600 C'tlFFri'44':Pl? - _ 

1602 F1-LEN-:P1:' : IF Fl=lTHEfJ 1440; IF i^9+r-i:-64THErJ 1630 


-161Q- 

_2T'R..:P81X 

1630 

OBleFSl . 1 

1610 


It' 40 

RE"'"UF'N 

•S^w 0 0 

DEFFN '24 

L'RTR 

SRVE Til 

2210 

L'LFFN'IS 

DRTR 

H T iii 

- 2220 

L'LFFN "S'E 

L'RTR 

SRVE T21 

2220 

DLFFN"20 

L'lTTR 

SRVE T41 


f >EEEfJ15-J 

j'3'40 

H=H I'l : on' 


H=iXF-ELLiJ4;— ■■ r .. E NEXT ...V EUfiEFF:" i^xasuaisso- 


_i;. .NE.VTT l..;FT:aNT H.i ^DflTfi SFiVE -ENE' 4 I^H-9THEM-23:60^L.riTR -SRV.E^r£N -"+J' :JjOT-- 
0 3'3'40 

5000 F'R I NT' 924 >: KS.- 1 .. E2'l ■' K.S.. 2 . E2’l •' KER 2 i' . B21 < K8. 4 B21 < KS.. 6 ? 

-6000 PRINT- Ell 2a>J^lXaC£L-2::x, B11C43£U-4.'X. Bll CX'ER-TO 
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